Strangeness dynamics and transverse pressure in relativistic nucleus-nucleus collisions 
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We investigate hadron production as well as transverse hadron spectra from proton-proton, 
proton-nucleus and nucleus-nucleus collisions from 2 A-GeV to 21.3 yl-TeV within two independent 
transport approaches (HSD and UrQMD) that are based on quark, diquark, string and hadronic 
degrees of freedom. The comparison to experimental data on transverse mass spectra from pp, pA 
and C-l-C (or Si-I-Si) reactions shows the reliability of the transport models for light systems. For 
central Au+Au (Pb-|-Pb) collisions at bombarding energies above ~ 5 A-GeV, furthermore, the 
measured transverse mass spectra have a larger inverse slope parameter than expected from 
the default calculations. We investigate various scenarios to explore their potential effects on the 
spectra. In particular the initial state Cronin effect is found to play a substantial role at top 
SPS and RHIC energies. However, the maximum in the /tt^ ratio at 20 to 30 A-GeV is missed 
by 40% and the approximately constant slope of the spectra at SPS energies is not reproduced 
either. Our systematic analysis suggests that the additional pressure - as expected from lattice QCD 
calculations at finite quark chemical potential fiq and temperature T- should be generated by strong 
interactions in the early pre-hadronic/partonic phase of central Au-|-Au (Pb-|-Pb) collisions. 
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I. INTRODUCTION 

The phase transition from partonic degrees of freedom 
(quarks and gluons) to interacting hadrons is a central 
topic of modern high-energy physics . In order to un- 
derstand the dynamics and relevant scales of this transi- 
tion laboratory experiments under controlled conditions 
are presently performed with ultra-relativistic nucleus- 
nucleus collisions. Hadronic spectra and relative hadron 
abundances from these experiments reflect important as- 
pects of the dynamics in the hot and dense zone formed 
in the early phase of the reaction. Here different quark 
flavors {u, d, s, c, 6) - as constituents of the final hadrons 
- provide valuable information on the early high den- 
sity phase due to the substantial different mass scales 
for the light (u, d), strange (s) and charm (c) or bottom 
(b) quarks. Furthermore, as has been proposed early by 
Rafelski and Miiller in particular the strangeness de- 
gree of freedom might play an important role in distin- 
guishing hadronic and partonic dynamics. 

Lattice QCD calculations at vanishing quark chemi- 
cal potential and finite temperature indicate critical en- 
ergy densities for the formation of a quark-gluon plasma 
(QGP) of ~ 0.7-1 GeV/fm^ 0, which might aheady 
be achieved at Alternating Gradient Synchrotron (AGS) 
energies of ~ 10 ^-GeV for central Au-|-Au collisions 



0,1^0]. Moreover, the recent lattice QGD calculations 
at finite quark chemical potential /i^ 7] show a rapid 
increase of the thermodynamic pressure P with temper- 
ature above the critical temperature Tc for a phase tran- 
sition (or cross-over) to the QGP. The problem is how to 
relate these numbers to observables from nucleus-nucleus 
collisions. Estimates based on the Bjorken formula 8] for 
the energy density achieved in central Au+Au collisions 
suggest that the critical energy density for the forma- 
tion of a QGP of ~ 0.7-1 GeV/fm'^ is by far exceeded 
during a few fm/c in the initial phase of Au-|-Au colli- 
sions at Relativistic Heavy Ion Collider (RHIC) energies 
. Accordingly, a phase transition (or cross-over) might 
already happen at much lower bombarding energies. 

Experimentally, relativistic nucleus-nucleus collisions 
have been studied at beam energies from 0.1 to 2 A-GeV 
at the Schwerlonen-Synchrotron (SIS), from 2 to 11.6 
A-GeY at the AGS and from 20 to 160 A-GeV at the 
Super Proton Synchrotron (SPS) jilic|- This program 
has been extended Q at RHIC with Au-|-Au collisions at 
invariant energies y/s from ^ 20 to 200 GeV (equivalent 
energies in a fixed target experiment: 0.2 to 21.3 A-TeV). 

Whereas many previous experimental studies have fo- 
cused on longitudinal rapidity distributions for hadrons, 
the transverse mass (or momentum) spectra of hadrons 
are presently in the center of interest. On the one hand 
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a significant suppression of high transverse momentum 
hadrons in Au+Au colhsions compared to pp is observed 
at RHIC energies of = 200 GeV [Hi El H El. 
This suppression is often attributed to the energy loss of 
hig hly energetic particles in a hot colored medium (QGP) 
[I^ lla|. In fact, the recent observation by the PHENIX 
[la, STAR O and BRAHMS [Tl collaborations, that 
a similar suppression is not observed in d+Au interac- 
tions at midrapidities at the same energy, supports this 
idea. On the other hand, the measured transverse mass 
{rriT = {Px + w^)^/^) spectra of hadrons (heavier than 
pions) 
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at AGS, SPS and RHIC energies show a substantial hard- 
ening in central Au+Au collisions relative to pp interac- 
tions (cf. [Tnl This hardening of the spectra (or 
increase of the inverse slope parameter T in J^l) for low 
transverse mass or momentum is commonly attributed to 
strong collective flow, which is absent in the respective 
pp or pA collisions. 

The authors of E3| have proposed to interpret the ap- 
proximately constant slopes above ~ 30 A-GeV - the 
'step' - as an indication for a phase transition. This inter- 
pretation is also based on a rather sharp maximum in the 
K+/TT+ ratio at - 20 to 30 A- GeV in central of Pb-|-Pb 
collisions (the 'horn' 18] ). However, the present body 
of experimental data and our theoretical understanding 
does not allow for a stringent conclusion based only on 
the 'step' and 'horn' phenomena. In this paper we will 
address some further theoretical issues, but would like to 
point out that more experimental data, for example on 
collisions in the energy range of 20 < ^/s < 200 GeV and 
on weakly interacting hadrons, will be needed to clarify 
these issues from the experimental side. 

We will demonstrate in this work, that neither the 
pressure needed to generate a large collective flow - to ex- 
plain the hard slopes of the spectra with a 'plateau' 
at SPS energies - nor the experimental sharp maximum 
in the K^/tt'^ ratio at ~ 20 to 30 A-GeV is produced 
in the present transport models by the interactions of 
hadrons in the expansion phase of the hadronic fireball. 
Some aspects of this work have been published in Refs. 
[l^ l20j before; however, we here extend the previous 
studies to RHIC energies and consider a wide variety of 
dynamical scenarios. In particular we will explore two 
distinct effects: i) initial state Cronin enhancement and 
ii) heavy resonance formation on the transverse mass 
spectra of pions and kaons at all bombarding energies 
from lower AGS to RHIC energies. 

In our studies we use two independent relativistic 
transport models that employ hadronic and st ring de- 
grees of freedom, i.e. UrQMD 0,53 and HSD [MEi. 
They take into account the formation and multiple 
rescattering of hadrons and thus dynamically describe 
the generation of pressure in the hadronic expansion 
phase. This involves also interactions of leading pre- 



hadrons that contain a valence quark (antiquark) from a 
primary 'hard' collision (cf. Refs. [i3)|23)- In the parton 
language, the two transport models employed do not in- 
clude gluon degrees of freedom explicitly nor gluon-gluon 
and giuon-quark (antiquark) interactions. We will use 
these models to test, whether the description of nucleus- 
nucleus reactions in terms of pre- hadronic, hadronic and 
string degrees of freedom and their interactions is suffi- 
cient to explain the experimental data. 

Our work is organized as follows: In Section II we 
briefly recall the ingredients of the HSD and UrQMD 
transport approaches and point out the extensions incor- 
porated in UrQMD (versions 2.0 and 2.1). In Section HI 
the calculated excitation functions for pions, mesons 
and hyperons for central Au-|-Au (or Pb-|-Pb) collisions 
from 2 to 21300 A-GeV are shown (in comparison to 
the experimental data) as well as selected particle ratios. 
Section IV is devoted to transverse mass spectra of pions 
and mesons from pp to pA and central AA reactions 
in the same energy range. In Section V we will discuss 
various alternative scenarios within the transport calcu- 
lations to explore their sensitivity to the hadron trans- 
verse mass spectra. Section VI is devoted to a study of 
the 'Cronin effect' in central Au-|-Au collisions, which we 
attribute to a 'pre-hadronic' initial state scattering effect 
that also shows up in pA reactions. Section VII includes a 
comparison of the thermodynamic Lagrange parameters 
T and /is extracted from the UrQMD transport model in 
the central overlap regime of Au-|-Au collisions with the 
experimental systematics on chemical freeze-out configu- 
rations in the T, /is plane as well as recent lattice QCD 
results. Section VIII closes this study with a summary 
and discussion of open problems. 



II. TRANSPORT MODELS - HSD AND URQMD 

We use two independent relativistic transport ap- 
proaches for our study that employ hadronic and string 
degrees of freedom, i.e. UrQMD (Ultra-relativistic Quan- 
tum Molecular Dynamics) '57, ^ and HSD (Hadron- 
String Dynamics) .23ii . .24 ^6] in order to obtain an es- 
timate of systematic uncertainties within the transport 
models. Both approaches take into account the forma- 
tion and decay of strings as well as the multiple rescat- 
tering of hadrons in the expansion phase. These models 
can be considered as appropriate descriptions at least for 
the final stage of an ultrarelativistic nucleus- nucleus col- 
lision, where the strongly interacting hadrons dominate. 
The open question is, however, if these transport mod- 
els also incorporate the proper degrees of freedom in the 
very early phase of heavy-ion collisions. We stress that 
transport calculations allow for a systematic study of the 
change in the dynamics from elementary baryon-baryon 
or meson-baryon collisions to proton-nucleus reactions 
or from peripheral to central nucleus-nucleus collisions 
in a unique way without change in parameters. This is 
of central importance since the assumptions of thermal 
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and chemical equilibrium - as used in hydrodynamical 
or statistical models - do not necessarily hold in all 
these reactions and the transport studies allow to explore 
the amount of (thermal or chemical) equilibrium reached 
Iks, J9, 3a]. 

The default UrQMD 1.3 transport approach [H ^ 
includes all baryonic resonances up to masses of 2 GeV 
as well as mesonic resonances up to 1.9 GeV as tabulated 
by the Particle Data Group [sj] . For hadronic continuum 
excitations a string model is used with meson formation 
times of the order of 1-2 fm/c depending on the momen- 
tum and energy of the created hadron. For the present 
study - and the particular aim to address also jet produc- 
tion and fragmentation at RHIC energies - we addition- 
ally have incorporated PYTHIA (v 6.1) HI] in UrQMD 
1.3 similar to HSD. This novel and extended version we 
denote as UrQMD 2.0. 

Moreover, to explore the physical consequences of 
higher mass states as implemented in the Relativistic- 
Quantum- Molecular-Dynamics model (RQMD) [s^, we 
have introduced in UrQMD additional high mass reso- 
nances that are explicitly produced and propagated in 
s- channel processes with invariant masses up to ^/s < 
3 GeV. We do not assume these states to be excited in 
baryon-baryon collisions. In practical terms: If the mass 
of the produced mesonic or baryonic resonance is above 
~ 2 GeV^, i.e in the continuum, the formed resonance 
state is not fragmented as a string object like in UrQMD 
1.3 or UrQMD 2.0. We recall that a string of invariant 
mass 2 < M < 3 GeV essentially decays to a few pi- 
ous, a vector meson -I- pions etc. Instead, this high mass 
state is treated as a quasi-particle that decays according 
to the branching ratios of the resonances with a mass of 
1.9 to 2.0 GeV (as implemented in UrQMD). We recaU 
that baryon resonances 1.9 GeV dominantly decay to 
the 7r + channel and with some lower probability to the 
p{uj) + N final state. In this way light meson emission 
is suppressed by about 25% compared to a string of the 
same invariant mass. We note, that this recipe is still in 
line with the available irp —^ irX data on pion multiplici- 
ties. Furthermore, since the high mass resonances decay 
isotropically in their rest frame, this leads to an enhance- 
ment of transverse flow (cf. Section V.E). Furthermore, 
we use isotropic angular distributions for elastic mcson- 
baryon interactions below Vs^^ = 6 GeV. In this work 
we will denote this scenario as UrQMD 2.1. 

In the HSD approach nucleons, A's, N*(1440), 
N*(1535), A, S and S* hyperons, S's, S*'s and Q's as 
well as their antiparticles are included on the baryonic 
side whereas the 0^ and 1~ octet states are included in 
the mesonic sector. Inelastic hadron-hadron collisions 
with energies above ~ 2.6 GeV are described by 
the FRITIOF model jHl (incorporating PYTHIA v 5.5 
with JETSET v 7.3 for the production and fragmenta- 
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tion of jets "S^) whereas low energy hadron-hadron colli- 
sions are modeled in line with experimental cross sections 

Both transport approaches reproduce the nucleon- 
nucleon, meson-nucleon and meson-meson cross section 
data in a wide kinematic range [l^ l2ll |2^. Therefore, 
we expect that the 'hadronic' scattering mechanisms are 
sufhciently under control. However, in nucleus- nucleus 
reactions further elastic and inelastic processes on the 
parton level might occur that - beyond the interactions 
of 'leading' quarks/diquarks - are not included in the 
studies below. We mention that (pQCD-like) partonic 
scattering is incorporated e.g. in the multi-phase trans- 
port model (AMPT) [s^l, which is currently employed 
from upper SPS to RHIC energies, but is not suitable for 
lower energies. 

III. EXCITATION FUNCTIONS 

Whereas the underlying concepts of the transport the- 
oretical models are very similar, the actual realizations 
differ considerably. A systematic analysis of results from 
both models and experimental data for central nucleus- 
nucleus collisions from 2 to 160 A- GeV has shown that 
the 'longitudinal' rapidity distributions of protons, pions, 
kaons, antikaons and hyperons are quite similar in both 
models and in reasonable agreement with available data 
[l9l |. We here extend the latter studies to RHIC ener- 
gies, but focus on total and midrapidity hadron yields as 
well as particle ratios. For completeness we include the 
results from Ref. |3 in Figs. Hfel 

A. Pion and strange particle yields 

Figure ^ shows the excitation function of 
TT+,Tr~,K+,K~ and A -I- S° yields from 5% cen- 
tral (AGS energies, SPS at 160 A-GeV and at RHIC 
energies), 7% central (20, 30, 40 and 80 A-GeV), 10% 
central for A -f 160 A-GeV Au-hAu (AGS and 

RHIC) or Pb-t-Pb (SPS) collisions in comparison to 
the experimental data from Refs. [s^ [s^. El] (AGS), 
HIM Hi 113 (SPS) and mill El (RHIC) for midra- 
pidity (left column) and rapidity integrated yields (right 
column). The solid lines with open squares show the 
results from HSD, the dashed lines with open triangles 
indicate the calculations with UrQMD 1.3 and the long 
dashed lines with crosses depict the UrQMD 2.1 results. 

We note in passing, that the corrections to the pion 
and kaon yields due to the weak decay of resonances have 
not been included in our calculations in line with the ex- 
perimental analysis of the NA49 Collaboration (cf. Refs. 
0,0). However, the experimental subtraction of week 
decays is problematic and might have an impact on the 
measured kaon to pion ratio as pointed out in Ref. |45j |. 
This uncertainty has to be kept in mind throughout the 
following analysis. 
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FIG. 1; (Color online) The excitation function of ■k'^,tt' -,K+,K- and A + E° yields from 5% central (AGS energies, SPS at 
160 A-GeV and at RHIC energies), 7% central (20, 30, 40 and 80 A-GeV), 10% central for A + S° at 160 A-GeV Au+Au(AGS 
and RHIC) or Pb+Pb (SPS) collisions in comparison to the experimental data from Refs. jsl IsHsl (AGS), [sl. IsgHiol ill 
(SPS) and |4^ l43l l4 ^ (RHIC) for midrapidity (left column) and rapidity integrated yields (right column). The solid lines 
with open squares show the results from HSD whereas the dashed lines with open triangles and the long dashed lines with 
crosses indicate the results from UrQMD 1.3 and UrQMD 2.1, respectively. The lower theoretical errorbars at RHIC energies 
correspond to the yields for 10% central events. 
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FIG. 2: (Color online) The excitation function of n^,-K and K mesons (as in FigCJ, but up to 80 A GeV on a linear 

scale. 



As already demonstrated in Ref. (19(] the agreement 
between the transport models HSD (solid lines) and 
UrQMD 1.3 (dashed lines) is rather good as well as with 
the data available. The exception are the pion yields at 
the highest AGS energy and lower SPS energies, which 
are overestimated by both approaches (cf. U3i)- This 
is demonstrated in more detail in Fig. |21 where the tt^ 
and K"^ yields (as in Fig. ^ are displayed in the energy 
range up to 80 A GeV on a linear scale. The overesti- 
mation of the pion yields (in in) in this energy range is 
slightly more pronounced for UrQMD 1.3 than for HSD 
especially at the higher SPS energy. On the other hand, 



both transport models fall slightly low in the K'^ yields 
at 30 and 40 A GeV. The results from the novel version 
UrQMD 2.1 (long dashed Hues) differ from UrQMD 1.3 
essentially in the following aspect: The pion multiplicity 
at AGS and SPS energies is reduced by about 25% (cf. 
Section II) due to the specific decay scheme of higher 
mass resonances (and is now in better agreement with 
the data). On the other hand, UrQMD 2.1 gives a lower 
amount of kaons and especially antikaons (in 47r) than 
the 'default' HSD and UrQMD 1.3 approaches. 

At RHIC energies the UrQMD 1.3 calculations show a 
~ 30% higher pion yield at midrapidity than HSD while 
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the differences between the Air integrated yields are of 
the order of 10%, both between the models and com- 
pared to the data. The implementation of higher mass 
states in UrQMD 2.1 (cf. Section II) - together with 
the jet production described by PYTHIA - leads to an 
overestimation of pions, kaons and hyperons at RHIC 
energies. 

Note, that the maximum in the A + yield at midra- 
pidity (lower left part of Fig. ^ for ~ 30 A-GeV is re- 
produced by both transport approaches. However, this 
maximum disappears in the An integrated yields in line 
with the data (cf. lower right part of Fig. ^J. Within 
the transport models the maximum in the A -I- E° yield 
at midrapidity can be attributed to a large extent to a 
change from baryon dominated dynamics at AGS ener- 
gies to meson dominated dynamics at the higher SPS 
energies. 



B. K+/tt+,K' /tt- and {A + T,°)/n ratios 

In Fig. 13 we present the excitation function of the 
particle ratios i^+/7r+, K^/tt' and (A-|- S°)/7r from 5% 
central (AGS energies, SPS at 160 A-GeV and at RHIC 
energies), 7% central (20, 30, 40 and 80 A-GeV), 10% cen- 
tral for A -f at 160 A-GeV Au-I-Au (AGS and RHIC) 
or Pb+Pb (SPS) collisions in comparison to the experi- 
mental data from Refs. 37, 38J (AGS), 9, 39, 41? ] 
(SPS) and m mm (RHIC) for midrapidity yields (left 
column) and rapidity integrated yields (right column). 
The solid lines with open squares show the results from 
HSD, the dashed hues with open triangles indicate the 
UrQMD 1.3 calculations whereas the long dashed lines 
with crosses correspond to the UrQMD 2.1 results. The 
deviations between the transport models and the data 
are most pronounced for the midrapidity ratios (left col- 
umn) since the ratios are very sensitive to actual rapidity 
spectra. The /tt^ ratio in UrQMD 1.3 shows a maxi- 
mum at ~ 8 A-GeV and then drops to a roughly constant 
ratio of 0.12 at top SPS and RHIC energies. UrQMD 2.1 
follows the results from UrQMD 1.3 at AGS and SPS en- 
ergies but then increases slightly up to RHIC energies. In 
case of HSD a continuously rising ratio with bombarding 
energy is found for the midrapidity ratios. This rise up 
to the full RHIC energy is essentially due to an underes- 
timation of the pion yield at midrapidity (cf. Fig. 1 in 
Ref. |48j). The Att ratio in HSD, which is approximately 
constant from top SPS to RHIC energies, is larger than 
the ratio from UrQMD 1.3 due to the lower amoimt of 
pion production and a slightly higher yield (cf. Fig. 

We note that the lower amount of pions in HSD is 
essentially due to an energy-density cut (in the local rest 
frame), which does not allow hadron formation when is 
energy density is above 1 GeV/fm'^. 

Despite this seemingly good description of the separate 
particle yields, the experimentally observed maximum in 
the i^r+/7r+ ratio at 20-30 A-GeV is not reproduced in 
the models. Due to strangeness conservation, the same 



arguments hold for the (A -I- S'')/7r ratio, where the pro- 
nounced experimental maxima are underestimated due to 
the excess of pions in the transport models at top AGS 
energies (for HSD) and above 5 A-GeV (for UrQMD 
1.3). Here UrQMD 2.1 performs better in the ratios due 
to the lower amount of pions produced by the high reso- 
nance decays. Since the yields are well described by 
all approaches (cf. Fig. ^ the deviations in the K~ /tt~ 
ratios at SPS and RHIC energies in UrQMD 1.3 indeed 
can be traced back to the excess of pions. This suggests 
that the entropy production, which is reflected in the 
pion to baryon ratio at these energies, is overestimated 
in UrQMD 1.3 (and HSD) above ~ 5 A-GeV in central 
Au+Au collisions. Note, that in UrQMD 2.1 the entropy 
production is reduced by the generation of high mass 
states in s-channels, which essentially decay by emitting 
a lower amount of pions. 

It has been discussed in the community that the 
maximum in the K'^/n'^ ratio might be due to isospin 
effects in Au-f Au (Pb-|-Pb) collisions arising from a dif- 
ferent rapidity dependence of u and d quarks in pp, pn 
and nn collisions. We here address this issue and show 
in Fig. 21 the influence of isospin effects on the kaon to 
pion ratio. The left part of Fig. 0] displays the excita- 
tion function of the (vr^) / (7r+) ratio from central Au-|-Au 
(or Pb-|-Pb) collisions calculated within HSD (solid line 
with open squares) and UrQMD 1.3 (solid line with open 
triangles). The right hand panel shows both the isospin- 
sensitive (K^) / {tt~^) ratio (open symbols) and the (ap- 
proximately) isospin-insensitivc 2{K^) / ((7r+) -I- (tt")) ra- 
tio (solid symbols). The HSD results are shown as 
squares connected by solid lines while the UrQMD 1.3 re- 
sults are represented as dashed lines with triangles. The 
larger solid symbols indicate the experimental results for 
{K+)/{tt+) from Refs. Hill El El and the large open 
triangles correspond to the ratio 2{K^) /{{tt^) + (vr^)) 
from Ref. E^]. As seen from Fig. 21 the calculated 
(as well as experimentally measured) isospin-insensitive 
ratio 2{K^) / {{tt^) + (tt^)) is only slightly lower than 
the isospin-sensitive (K^) / (tt^) ratio. This result is due 
to the approximately constant {■k~) / {■k'^) ratio above 
~ 6 A-GeV (l.h.s. of Fig. 2J) in the transport calcu- 
lations as well as in the data. We thus conclude that 
isospin effects cannot be responsible for the maximum in 
the {K+)/{-k)+ ratio (r.h.s. of Fig. Ej). 



IV. TRANSVERSE MASS SPECTRA 

Despite the relatively good agreement of data and 
transport calculations in the longitudinal direction [T^ . 
significant deviations arc found in the transverse direc- 
tion as pointed out in Ref. 20] before. This demon- 
strates that longitudinal rapidity distributions are less 
sensitive to the hadronic interaction channels employed 
in the transport models. We here present a more com- 
plete and systematic analysis up to the full RHIC energy 
of = 200 GeV and investigate alternative scenarios. 
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FIG. 3: (Color online) The excitation function of /■k'^ , K /tt and (A + Yf')/-K ratios from 5% central (AGS energies, SPS 
at 160 A-GeV and at RHIC energies), 7% central (20, 30, 40 and 80 A-GeV), 10% central for A + E° at 160 A-GeV Au+Au 
(AGS and RHI^or Pb+Pb (SPS) collisions in comparison to the experimental data from Refs. (AGS), 11 113,0 

(SPS) and |4^l43l l4^ (RHIC) for midrapidity (left column) and rapidity integrated yields (right column). The solid lines with 
open squares show the results from HSD whereas the dashed lines with open triangles and the long dashed lines with crosses 
indicate the results from UrQMD 1.3 and UrQMD 2.1, respectively. 
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FIG. 4: (Color online) The excitation function of the (7r~)/(7r"'") (l.h.s.) and {K'^) / {-k) ratios (r.h.s.) from central Au+Au 
(or Pb+Pb) collisions. Left part: the IfSD results are shown by the solid line with open squares, whereas the UrQMD 1.3 
calculations are displayed in terms of the solid line with open triangles. The experimental data have been taken from Refs. 

Umiilli^. Right part: the solid line with solid squares (triangles) (ci. F ig. corresponds to the HSD (UrQMD 1.3) 

calculations for (K^) / (tt^) (as well as the experimental data from Refs. |9ll3a. l4ll |42| ) while the solid line with open squares 
(tri angles) shows the calculated 2{K'^) /((tt^) + (tt^)) ratio in comparison to the experimental data (open triangles) from Ref. 



We start with benchmark tests for the transport calcula- 
tions for pp and pA reactions. 



A. pp and pA reactions 

Our results for the transverse mass (or momentum) 
spectra of pions and kaons from pp reactions are pre- 
sented in Fig. |31 The upper part includes the comparison 
of the differential cross sections for tt"*" , 7r~ and A'^ from 
pp reactions at 12 and 24 GeV/c from HSD (solid lines) 
and UrQMD 1.3 (dashed lines) with the data from Ref. 

The lower left part of Fig. [S] shows the transverse 
mass spectra at midrapidity for tt^ , and ^^(xO.l) 
from pp reactions at a bombarding energy of 160 GeV 
from HSD (soUd Unes) and UrQMD 1.3 (dashed Unes). 
The thin lines here correspond to fits with the experimen- 
tal slope parameters 163 ± 10 MeV for tt" , 172 ± 17 MeV 
for K+ and 164 ± 16 MeV for R- from Ref. jH. The 
lower right part of Fig. [31 displays the transverse mass 
spectra at midrapidity for 7r~ and K~ from pp reactions 
at = 200 GeV from HSD (soHd lines), UrQMD 1.3 
(dashed lines, multiplied by a factor 2), and UrQMD 2.0 
(2.1) (long dashed lines). The full symbols, furthermore, 
stand for the data from the STAR Collaboration |43 |. 
The description of the transverse mass spectra for pions 
and kaons is reasonably good in both approaches up to 
^/s = 200 GeV. We note in passing that the inclusion of 
PYTHIA (v 6.1) in UrQMD 2.0 and UrQMD 2.1 is es- 
sential in order to achieve a good description of the data 
since the contribution from jet fragmentation is essential 
at RHIC energies. Without jet production and fragmen- 



tation the UrQMD model (version 1.3) underestimates 
the high momentum spectra severely (similar to RQMD 

Mm)- 

We continue with pA reactions and show in the left 
part of Fig. the results for the inverse slope parame- 
ters T in (^) for tt^ , and mesons at midrapidity 
from pA reactions at 14.6 GeV/c (A=Be, Al, Cu, Au), 
whereas the right part includes the results at 450 GeV/c 
{A— Be, S, Pb). The inverse slope parameters from HSD 
are shown as open symbols connected by solid lines and 
from UrQMD 1.3 (2.0) by open symbols (with crosses) 
connected by dashed lines. The full symbols (l.h.s.) cor- 
respond to the midrapidity data {{yiab) — 1.5,1.7,1.9) 
from the E802 Collaboration ^ while the NA44 data 
^] at 2.4 < yiab < 3.5, pr < 0.84 GeV/c for K+,K- 
and at 2.4 < yiab < 3.0, pr = 0.3 1.2 GeV/c (full di- 
amonds) and 3.1 < yiab < 4.0, pr < 0.64 GeV/c (full 
squares) for tt"*" are given on the right-hand side. 

The overall reproduction of the transverse slope pa- 
rameters at AGS and SPS energies for targets from Be 
to Pb by both transport models is sufficiently accurate. 
The inverse slope parameters from UrQMD 1.3 and HSD 
differ for pions and antikaons by about 10% at AGS ener- 
gies. These differences can be traced back to the different 
production channels in the transport models at lower en- 
ergies: in HSD the mesons are essentially produced by 
string decays, which are strongly forward peaked in mo- 
mentum, whereas in UrQMD 1.3 they partly stem from 
high mass resonance decays with an isotropic decay in 
momentum space. At SPS energies the meson produc- 
tion in both models proceeds via string formation and 
decay and leads to very similar results. This also holds 
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FIG. 5: (Color online) Upper part: comparison of the differential cross sections for 7r"'",7r and Kg from pp reactions at 12 
and 24 GeV/c from HSD (solid lines) and UrQMD 1.3 (=UrQMD 2.0) (dashed lines) with the data from Ref. 'Ho']. Lower 
part (left): the transverse mass spectra at midrapidity for ■k~,K'^ and _ft'~(x0.1) from pp reactions at a bombarding energy 
160 GeV from HSD (solid lines) and UrQMD 1.3 (dashed lines). The thin lines correspond to fits of the experimental slope 
parameters 163 ± 10 MeV for vr", 172 ± 17 MeV for K+ and 164 ± 16 MeV for K' from Ref. @. Lower part (right): the 
transverse mass spectra at midrapidity for -k~ and K~ from pp reactions at ^/s = 200 GeV from HSD (solid lines) and UrQMD 
1.3 (dashed lines, multiplied by a factor 2) as well as UrQMD 2.0 (=UrQMD 2.1) (long dashed lines). The full symbols indicate 
the data from the STAR Collaboration ,4^. 



for RHIC energies where additionally jet production and 
fragmentation contributes with increasing bombarding 
energy. 

We note in passing that when including the initial state 
'Cronin effect' (cf. Section VI) the inverse slope param- 
eters T in Fig. increase by less than 5% even for the 
Pb-target at 450 GeV/c. At 14.6 GeV/c there is no effect 
(within statistics) for all targets from Be to Au. 



B. A+A reactions 

We continue with nucleus-nucleus collisions and recall 
that the experimental mx spectra (at midrapidity) from 
central C-l-C and Si-|-Si coUisions at 160 A-GeY (63| are 
well described by HSD and UrQMD 1.3. This holds for 
the 7r~, i4r+ as well as K~ spectra (cf. Fig. 1 in Ref. 



This situation changes for central collisions of heavy 
nuclei like Au+Au or Pb-|-Pb. In Fig. [7| we dis- 
play the calculated transverse mass spectra at midra- 
pidity from HSD (sohd Uncs) and UrQMD 2.0 (dashed 
lines) for 7t^,K'^ and K~{x0.1) from central Au-|-Au 
(Pb-fPb) reactions at 4, 6, 8, 11, 20 A-GeV (left part) 
and at 30, 40, 80, 160 A-GeV and ^/i = 200 GeV (right 
part) with the data (ful l and open symbols) from Refs. 
U3 In m in 2. At the lowest energy of 4 

A- GeV the agreement between the transport approaches 
and the data is still acceptable, however, severe devia- 
tions are visible in the spectra at top AGS, SPS and 
RHIC energies. Note that the tt^ spectra are reasonably 
described at all energies while the inverse slope T of the 



2 Note that all data from the NA49 Collaboration at 30 A-GeV 
have to be considered as 'preliminary' 
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FIG. 6: (Color online) Inverse slope parameters T for -k^jK^ and K~ at midrapidity from pA reactions at 14.6 GeV/c {A =Be, 
Al, Cu, Au) - left part and at 450 GeV/c {A =Be, S, Pb) - right part, from HSD (open symbols connected by solid lines) 
and UrQMD 1.3(2.0) (open symbols with crosses connected by dashed lines). The full symbols in the left part correspond 
to the midrapidity data {{yiab) = 1-5, 1.7, 1.9) from the E802 Collaboration Is^l . in the right part to the NA44 data [s^ 
at 2.4 < yiab < 3.5, pr < 0.84 GeV/c for K+,K- and at 2.4 < yiab < 3.0, pr = 0.3 ~ 1.2 GeV/c (full diamonds) and 
3.1 < yiab < 4.0, pt < 0.64 GeV/c (full squares) for 7r+. 



K transverse mass spectra is underestimated severely 
by about the same amount in both transport approaches 
(within statistics). The increase of the inverse slopes 
in heavy-ion collisions with respect to pp collisions, which 
is generated by rescatterings of produced hadrons in the 
transport models, is only small because the elastic meson- 
baryon scattering is strongly forward peaked (in HSD and 
UrQMD 2.0) and therefore gives little additional trans- 
verse momentum especially at midrapidity. 



V. TESTING ALTERNATIVE SCENARIOS 

In this Section we will investigate different 
hadron/string scenarios, which might contribute to 



the increase of the inverse slope parameter of the K 
mesons in central Au-fAu collisions and are not included 
in the 'default' versions of the transport approaches. 



A. In-medium potentials 

The question comes up whether the discrepancies 
shown in Fig. [7| might be due to conventional hadronic 
medium effects. In fact, the mx slopes of kaons and 
antikaons at SIS energies (1.5 to 2 A-GeV) were found 
to differ significantly [SJI- As argued in j2J| the differ- 
ent slopes could be traced back to repulsive potentials 
of kaons with nucleons, which lead to a hardening of the 
spectra, and attractive antikaon-nucleon potentials. 
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FIG. 7: (Color online) Comparison of transverse mass spectra at midrapidity from HSD (solid lines) and UrQMD 2.0 (dashed 
lines) for Tr^,K+ and K~{x0.1) from central Au+Au (Pb+Pb) reactions at 4, 6, 8, 11, 20 A-GeV (left part) and at 30, 40, 80, 
160 A-GeV and ^/s = 200 GeV (right part) with the data (full and open symbols) from Refs. |ia .36. .37.. .42 , .43 . .441 . 
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which lead to a softening of the K~ spectra. However, 
the effect of such potentials was calculated within HSD 
and was found to be of minor importance at AGS and 
SPS energies since the meson densities are compara- 
ble to or even larger than the baryon densities at AGS 
energies and above. 

Additional self energy contributions stem from in- 
teractions with mesons which become increasingly im- 
portant for higher bombarding energies; however, s-wave 
kaon-pion interactions are weak due to chiral symmetry 
arguments and p-wave interactions such as tt -I- if K* 
transitions are suppressed substantially by the approxi- 
mately 'thermal' pion spectrum. A recent study on the 
kaon potentials in hot pion matter gives kaon mass shifts 
of about —52 MeV and vector potentials of ^ -(-49 MeV 
[ssf for a pion gas at temperature T= 170 MeV. We have 
employed even slightly larger potentials in dynam- 
ical HSD calculations and achieved a hardening of the 
spectra by less than 10%. 



B. Overlapping strings 

Furthermore, we have pursued the idea of Refs. [sMIS^ 
that the K"^ spectra could be hardened by string-string 
interactions, which increase the effective string tension a 
and thus the probability to produce mesons at high rriT 
Is^ . A related suggestion has been put forward in 

In order to estimate the largest possible effect of string- 
string interactions we have assumed that for two overlap- 
ping strings the string tension a is increased by a factor 
of two, for three overlapping strings by a factor of three 
etc. Here the overlap of strings is defined geometrically 
assuming a transverse string radius i?s, which accord- 
ing to the studies in Ref. [53 should be Rs < 0.25 fm. 
Additionally an overlap of the strings in longitudinal di- 
rection is required which, however, does not involve any 
new parameter since the dynamics of the string ends is 
directly available from the HSD transport calculations. 
Based on these assumptions (and i?s=0.25 fm as in Refs. 

we find only a small increase of the inverse slope 
parameters at AGS energies, where the string densities 
are low (dashed lines in Fig. |S1 in comparison to the 
'default' calculations (sohd fines)). At SPS and RHIC 
energies the model gives a hardening of the spectra by 
about 15% which, however, is still too low in compari- 
son to the data in Fig. [7| Despite of the large number 
of 'initial' strings in central Au-|-Au collisions at SPS or 
RHIC energies, their average 'eigenvolume' (in the rest 
frame at break up) is low (~ 0.3 fm'^) compared to the 
eigenvolume of a nucleon (« 2 fm'^) due to the fast string 
fragmentation process. Consequently the critical density 
of color neutral strings may be about 7 times higher than 
the critical density of nucleons with respect to a phase 
transition induced by percolation [soj . 



T 1 1 1 1 f 

Au+Au, 11 A GeV, 5%, midrapidity 




T ' 1 ' 1 ' T 

Pb+Pb, 160 A GeV, 5%, midrapidity 




HSD, default (1) 

(l)+overl. strings 

(l)+parton elastic scattering 

(l)+isotropic decay of mB strings 

1 I I I I 



0.2 



0.4 





1 1 1 1 1 1 1 1 1 J 




Au+Au, s"'=200 GeV, 5 %, midrapidity 














K*0.1 




r , 1 , 


ri^i-'M 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

m^-m„ [GeV] 

FIG. 8: (Color online) The transverse mass spectra at midra- 
pidity for ■K~,K'^ and K~ (xO.l) from 5% central Au+Au 
(Pb-hPb) reactions at 11 A-GeV, 160 A-GeV and = 200 
GeV calculated for different scenarios in HSD: the solid lines 
correspond to the default calculations, the dashed lines to 
simulations with overlapping strings, the dotted lines reflect 
calculations including 'parton elastic scattering' and the dot- 
dashed lines arise from calculations with an isotropic decay 
of meson-baryon (mB) strings (see text). The experimental 
data are the same as in Fig. Q 



C. Parton elastic scattering 

We have examined another possibility which we de- 
note as 'parton elastic scattering'. To this end we have 
included elastic scattering of unformed hadrons in the 
initial phase of the reaction - before they actually evolve 
to their full hadronic size - employing the elastic pion- 
nucleon cross section per quark degree of freedom. The 
results of these simulations are shown in Fig. [S] by the 
dotted lines which give a much better reproduction of 
the experimental spectra at ^/s = 200 GeV, that are even 
slightly overestimated. The drawback on the other side is 
that the hadron rapidity spectra now are overestimated 
by ^ 20% relative to the default calculation (and the 
experimental data). The increase in the slope parame- 
ter especially for spectra at AGS and SPS energies. 



IS 




FIG. 9: (Color online) The rapidity distribution for 
■K~ ,K'^,K~ and protons from 5% central Pb+Pb reactions 
at 160 A-GeV. The solid lines correspond to the default HSD 
calculations, whereas the dashed lines show the calculations 
with an isotropic decay of strings from mB collisions (see 
text). The experimental data are taken from Refs. [9ll62|. 



however, is only moderate and a significant discrepancy 
remains with respect to the experimental spectra. 



D. Isotropic decay of strings from meson-baryon 
collisions 



Following a further idea of Rcf. we explore the 

possibility that the fragmentation of strings from mB 
collisions modifies in the dense and hot medium. An ex- 
treme modification amounts to assume an isotropic decay 
of mB strings - as employed in the RQMD model |33| 
- which leads to a reduction of longitudinal hadron mo- 
menta and an enhancement of their transverse momenta. 
The results of these simulations are displayed in Fig. |Slfor 
the transverse mass spectra of pions, kaons and antikaons 
(dot-dashed lines) . Apparently, this modified fragmenta- 
tion function leads to much better results in comparison 
to experiment (cf. dot-dashed lines in Fig. IS)). 

The drawback with this scenario, however, is that the 
approximate agreement with the longitudinal rapidity 
distributions for tt", , and protons is destroyed as 
demonstrated in Fig. I^for 5% central Pb-fPb reactions 
at 160 A-GeV. Here the solid lines correspond to the de- 
fault HSD calculations whereas the dashed lines stand for 
the calculations with an isotropic decay of mB strings. 
Obviously, the latter scenario leads to an enhanced stop- 
ping of the nuclei and larger meson production, which is 
not supported very much by the data. We admit, that 
a slightly more isotropic decay of the mB strings (with 
respect to the default calculations) would be favored by 
the experimental data; however, we do not aim at fitting 
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FIG. 10; (Color online) The transverse mass spectra at midra- 
pidity for -k' and if~(x0.1) from 5% central Au+Au 
(Pb-hPb) reactions at 11 A-CeV, 160 A-GeV and = 200 
GeV calculated for different scenarios in UrQMD: the dashed 
lines correspond to the default calculations (UrQMD 2.0), the 
solid lines to simulations with high mass baryon resonances 
(UrQMD 2.1). The experimental data are the same as in Fig. 
□ 



data here. We mention - without explicit representation 
- that the comparison shown in Fig. I^for central colli- 
sions at 160 A-GeV holds qualitatively also at AGS and 
RHIC energies. 



E. High mass baryon resonances (UrQMD 2.1) 

A further possibility to harden the meson and nucleon 
spectra in a coupled-channel transport approach is to in- 
troduce higher mass baryon resonances with masses > 2 
GeV as in RQMD 's^ , that decay isotropically and dom- 
inantly couple to the -nN channel and not to multipion + 
nucleon channels. As described in Section H, the version 
UrQMD 2.1 has been set up to incorporate this possi- 
bility and to explore its consequences. Since the 'mas- 
sive' resonances decay isotropically in their rest frame, 
this leads to enhanced transverse fiow as found in the 
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FIG. 11: (Color online) The rapidity distribution for 
■K~ ,K'^,K~' and protons from 5% central Pb+Pb reactions 
at 160 A-GeV. The dashed lines correspond to the default 
calculations from UrQMD 2.0, whereas the solid lines show 
the calculations with high mass baryon resonances (UrQMD 
2.1). The experimental data are taken from Refs. [gIIs^. 



previous Subsection. The penalty of enhanced hadron 
production (cf. Fig. is cured in UrQMD 2.1 since the 
high mass resonance decays produce 25% less particles 
as compared to a string of the same invariant mass. In 
this way one can produce larger hadron transverse mo- 
menta in the transport calculations without increasing 
the meson multiplicities on average. Such a 'prescrip- 
tion' fits the data in the AGS and SPS energy regime 
very well (similar to RQMD [ssLlsH ). However, the open 
question remains, if such 'states' could be identified in 
meson-baryon scattering and, furthermore, if they would 
survive in a dense hadronic medium. 

In Fig. ^1 we show the transverse mass spectra at 
midrapidity for 7r~,i4r+ and K^{'x.Q.\) from 5% central 
Au-FAu (Pb+Pb) reactions at 11 A-GeV, 160 A-GeV 
and ^/s = 200 GeV calculated for different scenarios 
in UrQMD: the dashed lines correspond to the default 
calculations (UrQMD 2.0), the solid lines to simulations 
with high mass baryon resonances (UrQMD 2.1). Indeed, 
the description of the pion and improves substan- 
tially in UrQMD 2.1 relative to UrQMD 2.0 (and HSD) 
at AGS and SPS energies. However, the penalty here is 
that now the spectra at the top RHIC energy are overes- 
timated at high transverse mass for pions and antikaons. 

Moreover, the improvement in the spectral slope for 
pions and kaons at AGS and SPS energies in UrQMD 
2.1 does not go along with an improvement in the rapid- 
ity spectra as shown in Fig. 1111 Though the tt^ rapidity 
distribution can be well described for 5% central Pb-|-Pb 
reactions at 160 A- GeV due to the reduction of pion pro- 
duction in the decay scheme of the high mass states, the 
rapidity spectra of mesons are also lower in UrQMD 



2.1 (solid lines) compared to UrQMD 2.0 (dashed lines) 
and the data. The proton stopping at midrapidity is not 
much affected by the introduction of the high mass states 
since they do not couple to the baryon-baryon channel. 

We conclude, that the recipe of high mass states - as 
incorporated in UrQMD similar to RQMD - does not 
solve the problem with the /tt^ ratio. 



VI. CRONIN ENHANCEMENT IN pA AND AA 
COLLISIONS 



The Cronin effect dates back to measurements of trans- 
verse hadron spectra from pA collisions at the SPS and 
ISR [6^ 1^ and relates to a (moderate) enhancement 
of transverse momentum spectra with increasing target 
mass A. Though the microscopic origin of this effect 
is not fully understood, most current interpretations are 
based on semi-hard initial state gluon radiation of the 
propagating quarks in the medium prior to a hard inter- 
action vertex (cf. e.g. |6,'t| ) . In this way the initial trans- 
verse momentum distribution of partons is extended to 
higher average (fc|^). We discard a further discussion of 
its physical origin and proceed with effective simulations 
following Ref. 25]. In the latter work it was suggested 
that such an enhancement of the intrinsic quark trans- 
verse momentum spread (k^) might be simulated in the 
transport approach by increasing the average transverse 
momentum spread of the quarks (fc^) with the number 
of previous collisions of primary nucleons Nprev as 



(4) 



(^t)pp(I + aNprev)- 



(2) 



The parameter a ~ 0.4 has been determined in Ref. |25j 
by comparing to the experimental data for d-|-Au colli- 
sions p^H^ at -^i = 200 GeV. This assumption is in line 
with an independent analysis within the Glauber model 
by Papp et al. 16^. Note, that Nprev is not a further pa- 
rameter of the model but calculated dynamically in the 
transport approach. 

Whereas the different alternative scenarios discussed in 
Section V relate to various model assumptions, which are 
not directly controlled by experimental data, the initial 
state Cronin enhancement is an established experimental 
fact. Consequently, it should be incorporated in 'future' 
transport approaches (as in Ref. jl^). For further details 
we refer the reader to the latter publication. 



A. A + A reactions with initial state Cronin 
enhancement 

Our results for the transverse mass spectra at midra- 
pidity for Tr~,K~^ and i^~(x0.1) from 5% central 
Au+Au (Pb-hPb) reactions at 11 A-GeV, 160 A-GeV and 
y/s — 200 GeV as calculated within HSD including the 
Cronin effect (0 (solid lines) are presented in Fig. [T21 
in comparison to the default calculations (dashed lines). 
The experimental data are the same as in Fig. [7| The 
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FIG. 12: (Color online) The transverse mass spectra at midra- 
pidity for ■K~,K'^ and K~{x0.1) from 5% central Au+Au 
(Pb+Pb) reactions at 11, 160 A-GeV and ^/s = 200 GeV 
calculated within HSD including the Cronin effect (solid 
lines) in comparison to the default calculations (dashed lines). 
The experimental data are the same as in Fig. |7| 



description of the spectra becomes rather good at RHIC 
energies, improves substantially at the SPS energy of 160 
A- GeV, but does not give any significant change at 11 
A- GeV. Consequently, the 'prehadronic' Cronin effect - as 
implemented via ||2Jl - is not responsible for the dramatic 
increase of the slopes at AGS energies. Our calcula- 
tions, however, show that a hardening of the transverse 
mass spectra should not only be attributed to the pres- 
ence of collective flow, but also to initial state semi-hard 
scattering of fast quarks in a dense medium; the latter 
effect becomes more important with increasing y^. A 
more detailed and quantitative analysis on the compe- 
tition of initial state Cronin enhancement and collective 
flow from pressure gradients at SPS and RHIC energies 
we postpone to an upcoming study. 



VII. INVERSE SLOPE PARAMETERS AND (pr) 

Our findings are summarized in Fig. Il^^l where the de- 
pendence of the inverse slope parameter T (defined by 
Eq. on -^i is shown and compared to experimen- 
tal dat a (p artly preliminary) from Refs. jlO, |3a Ell 
miM'STl for central Au+Au (Pb+Pb) coUisions (l.h.s.) 
and |4J,,6^,6^ for pp collisions (r.h.s.). The upper and 
lower solid lines (with open circles) on the l.h.s. in Fig. 
[TSl correspond to results from HSD calculations, where 
the upper and lower limits are due to fitting the slope T 
itself, an uncertainty in the repulsive K^-pion potential 
or the possible effect of string overlaps. The solid lines 
with stars correspond to the HSD calculations with the 
Cronin effect Q. The dashed lines with open upper tri- 
angles represent slope parameters from UrQMD 1.3, the 
dot-dashed lines with open lower triangles correspond to 
UrQMD 2.0 results, which are well within the limits ob- 
tained from the different HSD calculations without the 
Cronin enhancement. The dotted lines with crosses show 
the UrQMD 2.1 results that incorporate the high mass 
resonance states. 

The slope parameters from pp collisions (r.h.s. in Fig. 
I13|l are seen to increase smoothly with energy both in 
the experiment (full symbols) and in the HSD calcula- 
tions (full lines with open circles). The error bars re- 
lated to our calculations are due to an uncertainty in 
extracting the slope parameter from different ttit inter- 
vals. The UrQMD 1.3 results are shown as open upper 
triangles connected by the solid line and systematically 
lower than the slopes from HSD at all energies. When 
including jet production and fragmentation via PYTHIA 
in UrQMD 2.0 (dot-dashed lines with open lower trian- 
gles) the results become similar to HSD above -y/s = 10 
GeV demonstrating the importance of jets in pp reactions 
at high energy. 

Coming back to the slope parameters of mesons 
for central Au+Au/Pb+Pb collisions (l.h.s. of Fig. 1131) 
we find that the Cronin initial state enhancement indeed 
improves the description of the data at RHIC energies, 
however, does not give any sizeable enhancement at AGS 
energies. Here UrQMD 2.1 (dashed lines with crosses) 
with the high mass resonance states performs better for 
K'^ mesons, but overestimates the K~ slopes at AGS as 
well as upper SPS energies. 

An alternative way of representing the data is provided 
by the excitation function for the average transverse mo- 
mentum (pt) for pions (l.h.s.) and kaons (r.h.s.) in Fig. 
rPil for the different transport models in comparison to 
the data available. HSD: upper and lower solid lines 
with open circles, UrQMD 2.0: dot-dashed lines with 
open lower triangles and UrQMD 2.1: dotted lines with 
crosses; the data are from Refs. ^ 43, 44, 7(J E3- 
The dotted line with the open diamonds (r.h.s.) shows 
(pt) for K~ from UrQMD 2.1. Furthermore, the solid 
lines with stars correspond to the HSD calculations with 
the Cronin effect, which gives a good description of the 
average (px) in case of pions and significantly improves 
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FIG. 13: (Color online) Comparison of the inverse slope parameters T for K and K~ mesons from central Au+Au (Pb+Pb) 
collisions (l.h.s.) and pp reactions (r.h.s.) at midrapidity as a function of the invariant energy ^/s from HSD (upper and 
lower solid lines with open circles), UrQMD 1.3 (dashed lines with open upper triangles), UrQMD 2.0 (dot-dashed lines with 
open lower triangles), UrQMD 2.1 (dotted lines with crosses) with the data from Refs. .ICl, j36. ii. 42^ J3, J4j ^ for AA 
and 13, IgM |69| for pp collisions The upper and lower solid lines in the left diagrams result from different limits of the HSD 
calculations as discussed in the text while the solid lines with stars correspond to HSD calculations with the Cronin initial state 
enhancement J^J. 



the average (px) for kaons as a function of y/s (with re- 
spect to default HSD and UrQMD 2.0 calculations). The 
UrQMD 2.1 calculations overestimate {pt) for pions at 
all energies by 10 to 15 % as well as for K~ at AGS 
and SPS energies in line with the trends from Fig. 1131 
however, the average px for mesons is rather well 
described. 



VIII. THERMODYNAMICS IN THE T /is 
PLANE 

This still leaves us with the question of the origin of 
the rapid increase of the slopes with invariant en- 
ergy for central Au-|-Au collisions at AGS energies and 
the constant slope at SPS energies (the 'step'), which is 
missed in all transport approaches presently employed. 
We recall that higher transverse particle momenta either 
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FIG. 14: (Color online) The excitation function for the average values of pr at midrapidity for pions (left) and kaons (right) 
calculated within HSD (upper and lower solid lines with open circles) and UrQMD 2.0 (dot- dashed lines with open lower 
triangles), UrQMD 2.1 (dotted lines with crosses) in comparison to the data from Refs. 14^ liS. li^ l70l FH . The solid lines 
with stars correspond to the HSD calculations with the Cronin effect The dotted hne with the open diamonds (right part) 
gives (pt) for K' from UrQMD 2.1. 



arise from repulsive self energies - in mean- field dynamics 

- or from collisions, which reduce l ong itudinal momenta 
in favor of transverse momenta '72]. As shown above 
in Fig. ^1 conventional hadron self-energy effects and 
(pre-)hadronic and hadronic binary collisions are insuffi- 
cient to describe the dramatic increase of the slopes 
as a function of ^/s as well as the initial state Cronin en- 
hancement. Only when including high mass resonances 
in UrQMD 2.1 - with a particular decay pattern - we may 
end up with reasonable results for mesons, however, 
fail more severely for pions as well as antikaons. These 
failures of 10 to 15% indicate an additional mechanism 
for the generation of the pressure that is observed exper- 
imentally especially at AGS energies. 

We speculate - following the previous study (20j - that 
partonic degrees of freedom should be responsible for this 
effect already at ~ 5 A-GeV. Our arguments here are 
based on a comparison of the thermodynamic parame- 
ters T and /i^ extracted from the transport models in the 
central overlap regime of Au-f-Au collisions 28] with the 
experimental systematics on chemical freeze-out configu- 
rations in the T,^b plane. The solid line in Fig. [TBI 
characterizes the universal chemical freeze-out line from 
Cleymans et al. [tJ whereas the full dots with errorbars 
denote the 'experimental' chemical freeze-out parameters 

- determined from the thermal model fits to the experi- 
mental ratios - taken from Ref. |73l ]. The various sym- 
bols (in vertical sequence) represent temperatures T and 
chemical potentials /is extracted from UrQMD 1.3 trans- 



port calculations in central Au-f^Au (Pb+Pb) collisions 
at 21.3 A-TeV, 160, 40 and 11 A-GeV 28] as a function of 
the reaction time in the center-of-mass (from top to bot- 
tom). The open symbols denote nonequilibrium configu- 
rations and correspond to T parameters extracted from 
the transverse momentum distributions, whereas the full 
symbols denote configurations in approximate pressure 
equilibrium in longitudinal and transverse direction. 

During the nonequilibrium phase (open symbols) the 
transport calculations show much higher temperatures 
(or energy densities) than the 'experimental' chemi- 
cal freeze-out configurations at all bombarding energies 
(> 11 A- GeV). These numbers are also higher than the 
tri-critical endpoints extracted from lattice QCD calcula- 
tions by Karsch et al. lU] (large open circle) and Fodor 
and Katz (star with horizontal error bar). Though 
the QCD lattice calculations differ substantially in the 
value of fiB for the critical endpoint, the critical temper- 
ature Tc is closer to 160 MeV in both calculations, while 
the energy density is of the order of 1 GeV/fm^ or even 
below. Nevertheless, this diagram shows that at RHIC 
energies one encounters more likely a cross-over between 
the different phases when stepping down in temperature 
during the expansion phase of the hot fireball. 

This situation might change at lower SPS or AGS (as 
well as new GSI SIS-300 u^) energies, where for suffi- 
ciently large chemical potentials /i b th^ cross over should 
change to a first order transition (76j, i.e. beyond the tri- 
critical point in the {T,^b) plane. Nevertheless, Fig. E| 
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demonstrates that the transport calculation show tem- 
peratures (energy densities) well above the phase bound- 
ary (horizontal line with errorbars from Ref. 01 ) in the 
very early phase of the collisions, where hadronic inter- 
actions yield only a pressure p ~ 0.1 0.14e (with e 

denoting the energy density) . We recall that for massless 
partonic degrees of freedom we have p ^ e/3 and conse- 
quently a much higher pressure. This argument is well 
in line with the studies on elliptic flow at RHIC energies, 
that is underestimated by ~ 30% at midrapidity in the 
HSD approach and by a factor of ~ 2 in the UrQMD 
1.3 model f7% . Only strong pre-hadronic/partonic inter- 
actions might cure this problem. 

As shown in Ref. flS\, in order to describe the elliptic 
flow seen experimentally at RHIC in a parton cascade 
model, one has to employ parton cross sections up to 45 
mb. However, such strong interactions are incompatible 
with perturbative QCD, which gives cross sections that 
are lower by more than an order of magnitude |7l|. We 
speculate that such strong nonperturbative interactions 
on the partonic level are responsible for the large pressure 
generation in the very early phase of intermediate energy 
nucleus- nucleus collisions. 



IX. SUMMARY 

Summarizing this work, we point out that baryon 
stopping and hadron production in central Au-|-Au (or 
Pb-|-Pb) collisions is quite well described in the indepen- 
dent transport approaches HSD and UrQMD 2.0. Also 
the 'longitudinal' rapidity distributions of protons, pi- 
ous, kaons, antikaons and hyperons are similar in both 
models and in reasonable agreement with available data. 
The exception are the pion rapidity spectra at the high- 
est AGS energy and lower SPS energies, which are over- 
estimated by both models J3}. As a consequence the 
HSD and UrQMD 2.0 transport approaches underesti- 
mate the experimental maximum of the K'^ /tt'^ ratio at 
~ 20 to 30 A-GeV and UrQMD underestimates also the 
data at RHIC energies. This also holds for the new ver- 
sion UrQMD 2.1 introduced here additionally, which in- 
corporates high mass resonance states above 2 GeV with 
a particular decay pattern. This recipe cures the problem 
of the overestimated pion yield at AGS and SPS energies 
but also reduces the yield and therefore has only 
little effect on the if +/7r+ ratio. 

We have found that the inverse slope parameters T 
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for mesons from the HSD and UrQMD 2.0 transport 
models are practically independent of system size from 
pp up to central Pb+Pb collisions and show only a slight 
increase with collision energy. The calculated transverse 
mass spectra are in reasonable agreement with the ex- 
perimental results for pp reactions at all bombarding en- 
ergies investigated as well as central collisions of light 
nuclei (C-l-C and Si-hSi) (cf. Ref. ^20]). The rapid in- 
crease of the inverse slope parameters of kaons for central 
collisions of heavy nuclei (Au-|-Au or Pb-|-Pb) found ex- 
perimentally in the AGS energy range, however, is not 
reproduced by both models in their default version (see 
Fig. ESI). 

We have, furthermore, discussed a couple of scenarios 
in Section V to improve the description of the experi- 
mental data. However, none of the alternative scenarios 
leads to a fully satisfactory reproduction of all the exper- 
imental data. It is worth pointing out that a substantial 
hardening of the slope results from an initial state 
Cronin enhancement at higher SPS and RHIC energies, 
that should not be attributed to collective flow (or initial 
pressure). This Cronin enhancement ~ as known experi- 
mentally from p + A reactions at SPS energies - does not 
signal a phase transition to a quark-gluon-plasma. 

We thus have argued - in comparison to lattice QCD 
calculations at finite temperature and baryon chemical 
potential fiB from Refs. |7| and i7j| as well as the exper- 
imental systematics in the chemical freeze-out parame- 
ters (cf. Fig. 115(1 - that the missing pressure should be 
generated in the early phase of the collision by nonpertur- 
bative partonic interactions because the strong hadronic 



interactions don't produce it in the later stages. 

The interesting finding of our analysis is, that pre- 
hadronic/partonic degrees of freedom seem to play a sub- 
stantial role in central Au-|-Au collisions already at AGS 
energies above ~ 5 A-GeV. This finding is in line with the 
'observation' that the excitation function of the proton 
elliptic flow V2 also indicates a softening of the equation of 
state with increasing energy at 4-5 A-GeV in Au-|-Au col- 
lisions [t^Is^,!^- However, further detailed experimen- 
tal and theoretical studies also on collective flow prop- 
erties from nucleus-nucleus collisions will be necessary 
to unravel the various aspects of rclativistic many-body 
systems. 
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